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Abstract 

The effects of uniaxial strain on the structural, orbital, optical, and magnetic 
properties of LaMnOa are calculated using a general elastic energy expres- 
sion, along with a tight-binding parameterization of the band theory. Tensile 
uniaxial strain of the order of 2 % (i.e., of the order of magnitude of those 
induced in thin films by lattice mismatch with substrates) is found to lead 
to changes in the magnetic ground state, leading to dramatic changes in the 
band structure and optical conductivity spectrum. The magnetostriction ef- 
fect associated with the Neel transition of bulk(unstrained) LaMnOs is also 
determined. Due to the Jahn- Teller coupling, the uniform tetragonal distor- 
tion mode is softer in LaMnOs than in doped cubic manganates. Reasons 
why the observed (vrvrO) orbital ordering is favored over a (vrTrvr) periodicity 
are discussed. 
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I. INTRODUCTION 



The 'colossal' magnetoresistive manganese perovskites have been a focus of recent 
attention.S Since most of the technological applications require thin films on substrates, 
it is important to understand the effects of strains induced by substrates. Because Mn Cg 
electrons, which determine important physical properties of these materials, are coupled to 
the lattice degrees of freedom through the Jahn- Teller (JT) coupling, it is expected that 
uniaxial or biaxial strains are important and that even relatively small strains may result 
in observable effects on the electronic properties of these materials. Recently, the effects 
of substrate-induced strains on the properties of thin manganate films have been studied 
experimentally.lili It is indeed found that the ferromagnetic and metal-insulator transition 
temperature, Tc, depends sensitively on the volume-preserving uniaxial strains, as do the 
magnetic anisotropy, magnetoresistance, and charge ordering transition. 

In this work, we study the effects of uniaxial strains in LaMnOa, which is the parent 
compound of the doped manganese perovskites. Our motivations are to further clarify 
the properties of this interesting compound and to test techniques and fix parameters so 
the more complicated behaviors of the doped compounds may be studied. At very high 
temperatures, bulk LaMnOs exists in cubic perovskite structure, but at T < 750 K, it has a 
static (tt, tt, 0) 3a;^ — r^/Sy^ — type Jahn- Teller distortion. @ It also has a uniform tetragonal 
distortion, which originates from the coupling of the staggered and uniform distortions by 
an anharmonic elastic energy.! This coupling implies that the substrate-induced strain may 
affect the orbital ordering. 

In Ref. following Kanamori, a model with harmonic Mn-0 and Mn-Mn elastic forces 
and a local anharmonic energy term was used to study lattice distortions in this material. 
It predicted {it, it) type ordering in xy plane as observed in LaMnOa. However, according 
to this model, (vr, vr, 0) and (tt, tt, tt) type orderings have the same energy. The reason why 
(tt, tt, 0) ordering is favored has not been understood so far. Bulk LaMnOa has so called A- 
type antiferromagnetic (AF) ordering below 140 K, in which spins align parallel in xy plane 
and antiparallel along z direction.§ This peculiar spin ordering pattern is the result of the 
exchange interaction between Mn ions which depends on the Cg orbital ordering as explained 
in Ref. ^ Indications of a coupling between orbital ordering and magnetic ordering have 
also been found in a recent X-ray resonant scattering experiment.il 

In this paper, we present more general expressions for the elastic energy than those 
in Ref. |[ We use these to study the ground state energy and distortions in bulk state, 
and the effects of uniaxial strains in thin films. We also present a simple model of the 
magnetic interaction depending explicitly on the orbital states, and we use this to study 
the magnetostriction effects in bulk state and the change of the magnetic interaction due 
to the strains in thin films. We examine the changes in the band structure and optical 
conductivity due to the strains, using a nearest-neighbor tight binding approximation. We 
also compare the energies of (vr, vr, 0) and (vr, vr, vr) type orderings, and examine why the 
observed (vr, vr, 0) distortion is favored over the (vr, vr, tt) type distortion. Our calculations 
suggest that 2 % tensile strain can change the observed A-type (layered) antiferromagnetic 
state into a purely antiferromagnetic state. This change would induce large changes in band 
structure and optical conductivity spectrum, which we calculate. Finally, we show that the 
magnetostriction effect at the Neel transition is large. 
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The rest of the paper is organized in the following way. Section II introduces a model 
of elastic energy, Sect. Ill a model of the magnetic interaction and magnetostriction effect, 
and Sect. IV a tight binding model for band structure and optical conductivity. Section 
V compares the ground state energies of (vr, 7r,0) and (7r,7r,7r) type distortions. Section VI 
presents the results. Section VII summarizes the main conclusions. In Appendices A and B, 
we show how we determine general expressions of the elastic energy, and how we determine 
the parameters, respectively. 



To study the effects of uniaxial strains in thin films, we first need to understand the strains 
already present in bulk state. The elastic energy depends on three important variables: the 



consider are defined in the following way: In the ideal cubic perovskite structure with lattice 
constant Oq, the Mn ions are located at Qqi, and oxygen ions at ao(i + a/2), where i^, iy, and 
iz are integers and a=x, y, and z. We write the displacement of Mn at aoi as ao{e-i + 5^), 
where 5^ represents nonzero-wavevector Mn-ion displacement, e = CxxX + Cyyy + CzzZ, and Cij 
is the conventional strain tensor referred to the ideal cubic perovskite lattice (we need only 
the diagonal components). The displacement of oxygen at ao{i + a/2) along Mn-Mn axis 
is ao[e ■ (i + d/2) + u?a], where represents 0-ion displacement with nonzero wavevector. 

Figure 1 shows 6^ and u^'^'^. We assume that we have already minimized the elastic energy 
over the displacements of O ions perpendicular to the Mn-Mn axis and the displacement 
of La ions. Therefore, even though we do not consider the buckling of the Mn-O-Mn bond 
explicitly, its effect is implicitly included in our harmonic and anharmonic elastic constants 
below. 

We treat the elastic energy due to these strains in up to cubic anharmonic terms, and the 
coupling of the strain to the Mn Cg orbital state by the Jahn- Teller coupling. In this paper, 
instead of representing the elastic energy in terms of spring constants between different 
ions as done in Refs. ^ and |^, we will introduce a more general and, we hope, more useful 
formulation in which long wavelength lattice distortions are treated via conventional elastic 
theory, while the short wavelength atomic motions are treated explicitly. 



The energy per Mn ion due to uniform strains, e^x, ^yy, and e^z, is most conveniently 
written in terms of the following combinations: 



II. MODEL OF LATTICE ENERGY 



A. Overview 




B. General energy expressions for (7r,7r, 0) and (7r,7r, vr) distortions 



Q2U = ^(C; 




(1) 



(2) 
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We have 



Qsu ^j-^^(^zz (^xx ^yy)' (3) 



#^ = ^i^BQL + ^^*(QL + Q3J, (4) 



where the bulk modulus Kb = (3.o(cii + 2ci2) = Saoc^, the Jahn- Teller shear modulus 
K* = ao(cii — C12) = 2aoc*, and Qj are the usual elastic constants. 

In addition to these uniform strains, we consider staggered distortions with wave vector 
k = (tt, 71, 0) or (tt, tt, tt). We represent the amplitudes of the distortions by aoSj: for Mn ions 
and ooM^'^'^ for O ions. Translational symmetry implies that the uniform distortion and 
the staggered distortion do not couple with each other up to the second order. Therefore, 
the harmonic elastic energy due to the staggered distortions simply adds to the harmonic 
uniform strain energy. In Appendix A, we use symmetry arguments to obtain the general 
forms of the elastic energies due to 5^ and u^'^'^ for (vr, vr, 0) and (vr, tt, tt) distortions. In the 
Jahn- Teller coupling energy Ejt, which we will introduce later, only the lattice distortions 
which have even parity about Mn ions appear. Therefore, Ejt does not depend on the Mn- 
distortions ^^'^'^ for either k = (vr, vr, 0) or A; = (tt, tt, tt), or the z-direction oxygen distortion 

m| for k = {it, tt, 0). The energy cost of the relevant distortions is most conveniently written 
in terms of 

Q2s = "~|(^^^ ~ ^s?/)' (^) 

where Vsa is the (ttttO) or (tttttt) amplitude of = — The energy of a staggered 

distortion. Eg, depends upon the ordering wavevector and, for the distortions we consider, 
is 

E,(tT,TT,0) 1.^ ^2 ^2 



7;K2sQL + 7;K,sQL, (7) 



Nmu 2 --^'^ 2 

%^ = k(QL + QL), (8) 

where the K2S, K^g, and Kg are the elastic constants defined in Appendix A, and arise mainly 
from the Mn-0 bond stretching mode.i 

We compare the sizes of K*, K2S, K^g, and Kg. This comparison is important to deter- 
mine the ground state distortions. A uniform strain changes more bonds than a staggered 
distortion. For example, within a harmonic nearest-neighbor approximation, Q2S and Q^g 
modes involve only the Mn-0 bond, but Q2U and Qsu modes involve both Mn-0 and Mn-Mn 
bonds. Therefore, uniform modes have larger elastic moduli than staggered modes, which 
remains true when reasonable further neighbor interactions are included. Therefore, we ex- 
pect K* > K^g, K2g, Kg. Our analysis of the general expression of the elastic energy given 
in Appendix A shows K^g > K2g. 
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The staggered lattice distortion is caused by the Jahn- Teller coupling to the Mn Cg orbital 
state. The Mn Cg electron state on site a^i is represented by 



|%>= cosOMz^ -r^ > +sin%|a;^ - >, 



(9) 



where \3z'^ — > and \x'^ — y'^ > are the two linearly independent Cg orbitals on the site. 
The cases of interest here are two-sublattice distortions. We represent the orbital states of 
Cg electrons on these two sublattices as 6i and 62- Then the JT energy iJl'i 



"9 

E 



JT 



N, 



Mn 



[C0S2^1(Q3« + Q3s) + sin2^i(Q2« + Q2s) + cos 2^2 (Qs^ - Q3s) + sin2^2(Q2« - Q2s)] 



2 2 



(10) 



To represent the energy only in terms of lattice distortions, we minimize the above Jahn- 
Teller energy with respect to the orbital states 9i and 62. We obtain 



cos 2^5^^ 
sin2^;"^" 



3s 



{Q2u±Q2sy+{Q3u±Q3s)''' 

Q2U ± Q2S 
{Q2u±Q2sy+{Q3u±Q3sy'' 



(11) 
(12) 



and 

EjT 



Mn 



■In 



-A 



iQ2u + Q2sy + {Q3u + Q3sY + V(Q2« " Q2sY + (Q3« - Q3sY 



(13) 



We also include the largest anharmonic energy, which is the one between the nearest 
neighbor Mn-0 pair. It is given by 



Y 



(14) 



Total energy is the sum of the terms considered so far: 



rpelastic I? i I? i Z? i Z? 

^tot — I^u + J^s + ^JT + ^anh- 



(15) 



By minimizing Efj^l^^^'^, we find the distortions induced by the JT coupling and anharmonic 
energy terms, which will be discussed below for k = (tt, it, 0) and k = (tt, tt, it) distortions 
separately. 



C. Energy minimization for (vr, vr, 0) distortion 

We minimize Eu + Es + Ejt and then treat Eanh as a perturbation, since we expect and 
will show below that the anharmonic term is small compared to the harmonic terms. We 
find that in the ground state of E^ + Eg + Ejt, the distortion mode which has the smallest 
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modulus among Q2U, Qsu, Q2S, and Q^s is non-zero, and all the other distortion modes 
vanish. Since, we found K* > K^g > K2S, the ground state of + + Ejt is 

Qss = Q2U = Qzu = 0, (16) 

= "4:^' ^^^^ 

^1,^2= vr/4,37r/4. (19) 

It is noteworthy that in this ground state the Mn lattice itself preserves cubic symmetry, 
and only oxygen ions make staggered distortions. 

We next study how the anharmonic energy term changes the above ground state. We 
represent Eanh in Eq- (0) by 5^, m^'^'^, and exx,yy,zz- Direct expansion (or symmetry argu- 
ment) shows that each of 5g and m| appears only as the second order in Eanh, which implies 
that these distortions remain zero unless the coupled uniform strains exceed certain values 
and the lattice becomes unstable. After representing EanhiT^T^O) in terms of Qiu,2u,3u,2s,^ 
the same argument shows Qss = Q2U = 0. By taking the largest remaining term, we obtainll3 

%^-AQL(Qi.-^Q3.). (20) 

The total energy for k = {it, vr, 0) distortion which we will minimize to find strains in 
bulk state or in thin films is 

^tot(vr7r0) = {Eu + E, + Ejt + ^anh)/A^Mn 

= ^KbQI + Ik*QI + \k,sQI - ^X^Ql + QL + AQl (^Q^^ - -^Qa^) • (21) 

In bulk state, there is no external constraint. When we minimize £'tot(7i"7i"0) in a leading 
order in A to find the lattice distortions in the bulk state, we obtain 

for 

03, = ft- = 0, (25) 
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The results show that the observed uniform tetragonal distortion is due to the anharmonic 
term which couples staggered and uniform distortions. LaMnOa expands upon heating, 
which implies A < 0. Therefore, above result indicates Qsu < 0, which is consistent with 
the observed distortion in bulk LaMnOs. Since K* — is order of magnitude smaller than 
K2S1 Qsu is order of magnitude larger than Qin- Up to order of A, the anharmonic term 
does not change Q2s- 

We note "Qsu mode softening" : When we have lattice distortion Q2S = ^3/2-^, the 
JT coupling term in Eq. (pT]), expanded about small Qsu, effectively reduces Qsu niode 
modulus by —K2s- Therefore, when the anharmonic term induces Qzui the restoring spring 
constant is K* — K2S rather than K* . This is the reason why we have K* — K2S in the 
denominator of Eq. (|26|) and have a relatively large Qzu{^ 35 % of Q2S from crystallography 
data). Thus we expect the shear modulus corresponding to the Q^u distortion to be much 
smaller in LaMnOs than in doped cubic manganates. 



D. Energy minimization for (7r,7r, vr) distortion 



By applying similar considerations to (vr, vr, vr) distortion, and using the condition K* > 
Ks, we find that the degenerate ground states of Eu + Es{ tttttt) + Ejxi^TTTTTT^ are 



Q2U — Qsu — 0, 



/3 A 

Q2S = 1/ sin26'i. 



Q 



3s 



E --3^ 
^2 = vr-0i, 



COS 2^1, 



(27) 
(28) 

(29) 

(30) 
(31) 



where di is an arbitrary angle between and tt. After we include the same anharmonic 
energy and apply the same arguments used for the (tt, tt, 0) distortion, we find the total 
energy expression for the (vr, vr, vr) distortion which we will minimize further is 



Etot (tttttt) = [Eu + ^.(tttttt) + (tttttt) + Ean/, (tttttt)] /A^, 



Mn 



■is) 



(32) 
(33) 



2 2 



{Q2u + Q 



2s 



(Qiu + Qss)^ + JiQ2u - Q2sf + {Qiu - Qisf 



Q'isiQiu — -^Qd,u) + Q\s{Qiu + -^Qiu) — V2Q2sQ3sQ 



+A 

When we minimize Etot{T^T^T^) ■, we obtain 

j^MIN (^^\ 
J^tot ( TTTTTT) 



3 A^ 



(34) 
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for 



Q..- -^£1 + 0(^3), (35) 

l^^+OiA% (36) 

^ 0, (37) 

^^'^ +0{A'). (38) 




The other two physically equivalent distortions obtained by permuting x, y, and z from 
above results are also degenerate ground states. The ground state energies for the (vrvrO) 
and (tttttt) distortions will be compared in Sect. V. 



E. Energy minimization for strained films 

In this paper, by uniaxial strain, we mean a tetragonal strain with axis perpendicular to 
film plane. This strain can be applied by growing epitaxial films on square lattice substrates 
with lattice parameters different from the xy plane lattice parameter for bulk LaMnOs. For 
thin films with uniaxial strains, we calculate the lattice distortions in the following way: 
We assume that (ttttO) distortion pattern is favored even in strained films. We also assume 
perfect epitaxy; therefore Cxx and Cyy are determined by substrates. Then Q2s and Czz for 
given strains can be found by minimizing Etot about these distortions, which also give Qiu 
and Qsu. From Q2s and Q^u, we find the Cg orbital states, 9i and 02, by 

cos(2^^i) = ^ , (39) 

V Qh + Qlu 

62 = 71-61. (40) 

Since the external strain changes Q2S and Qsu, it also changes the JT splitting, given by 

6AEjT = VQX^QIs + Qlu - V^\\l{Q2sY + {QZY- Parameters of the model, i.e., Kb, 
K*, K2S1 A, and A, are determined from experiments, as explained in Appendix B. The 
determined parameter values are shown in Table I. 



III. MODEL OF MAGNETIC INTERACTION AND MAGNETOSTRICTION 

EFFECTS 

A. Magnetic interaction 

In this section, we present a model describing magnetism in bulk and strained films 
of LaMnOs. The superexchange interaction between Mn ions in LaMnOa depends on the 
electron orbital overlap, particularly the overlap between Mn Cg orbitals and O p orbitals.liil 
It is also argued that t2g electrons also contribute antiferromagnetic interaction. We build 
a simple model below, which incorporates these two contributions of exchange interactions. 
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We calculate the superexchange interaction due to the Cg electrons using a similar method 
as in Ref. ^ We assume Cg spin is always parallel to t2g spin at each site due to the strong 
Hund's coupling. For the two Mn ions, one at i and the other at i + z, we find 

Je,{Oi, 62) = -Jf (sin^ 6^ cos^ 62 + cos^ Oi sin^ 62) , (41) 

where Jp is a positive parameter of the model. This is, in fact, equivalent to the model in 
Ref. 1^, if we assume that the state with two holes on the intermediate oxygen ion (which 
was considered in Ref. ^ requires an infinite energy. Equation (|4l|) shows that when one 
of 6*1 and 6*2 is zero (i.e., |3z^ — > state) and the other is tt/2 (i.e., — y"^ > state), 
Jeg(6'i,^2) is most ferromagnetic due to the maximized hopping between the filled orbital 
on one site and the empty orbital on the other site. When 61 = 62 = 01 9i = 62 = 0, 
Je (^i! ^2) is least ferromagnetic, since the hopping between the filled and empty orbitals on 
the two sites vanishes. The t2g superexchange is expected to be independent of 61 and 6*2, 
and always antiferromagnetic, so we set Jt2g = Jaf- Therefore, along the z direction the 
total exchange interaction is 

Jz = Jt2g + Jeg 

= Jaf - Jf (sin^ 61 cos^ 62 + cos^ 61 sin^ 62) ■ (42) 

The sign of the total superexchange is determined by the competition between the Cg fer- 
romagnetism and t2g antiferromagnetism. Along x and y directions, proper rotations result 
in 

J. = Jaf - Jf {sin^iO, - ^) cos2(^2 - f ) + cos\9^ - ^) sin2(02 - f )) , (43) 

/ 2tt 2tt 2tt 2tt \ 

Jy = Jaf - Jf (sin^e, + —) cos2(^2 + y) + cos\e, + — ) sin2(^2 + y)) • (44) 

Using the condition 62 = n — 9i obtained before, we find 

J^ = Jy = J^y = Jly + J* cos 4^1, (45) 
= J° + J*cos4^i, (46) 

where J^y = Jaf — 5Jf/S, J° = Jaf — JfI^i and J* = Jf/4:. It explicitly shows how the 
magnetic coupling depends on the orbital states. 



B. Magnetostriction effects 

In this section, we present our model of magnetostriction effect in bulk LaMnOs. We 
expect magnetostriction effect for the following reason: Since J^, Jy, and J^ depend on the 
orbital states di and 62, once certain magnetic ordering occurs, the magnetic ordering, in 
turn, will change 61 and 62 to gain further magnetic energy. Through the JT coupling, this 
change in the orbital states can cause the change in the JT strains. Magnetism and lattice 
distortion are coupled through the Mn Cg orbital degree of freedom. 

To estimate the size of the magnetostriction effects, we add to Eq. (pll) the term 
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|^ = -4(2|J.,(^i)l + l^.(^i)l), (47) 

which represents the mean-field T=0 magnetic energy, and find the changes in the orbital 
states. For A-type antiferromagnetic state, after dropping constant terms, we obtain 

^ = 4J*cos4^i. (48) 

By adding Emag/NMn to the total energy, we can find the extra structure at T=0 due to the 
magnetic order. 

Using the Landau free energy method, we examine whether the magnetostriction effect 
makes the phase transition first order or not. From the energy gain due to the magnetic 
order at T=0, and the mean field estimate of T^, we obtain the Landau free energy per site, 
/: 

/(^i, Q2s, Q3«, m) = 2[T - T,{ei)]m^ + T,{ei)m^ + E,iasUc{Oi, Qm, Q2s, Qsu), (49) 

where m = M/M(T = 0) is the normalized magnetization and EeiasUci^i, Qiu, Q2s, Qsu) is 
the total elastic energy obtained in Sect. II. Magnetic ordering temperature, Tc{9i), can be 
obtained by molecular field theory:lll 

TM = 2{A\J,y{d,)\ + 2\Ue,)\). (50) 

We can find the order of the phase transition and the size of the magnetostriction effect in 
the following approximation: Since Qiu is not directly coupled to 9i, we expect the change 
in Qiu is small, which can be seen from the numerical results in Table II. Therefore, we 
neglect Qiu dependence and obtain Eeiastic{Oi, Q2s, Qsu) , which we minimize further about 
Q2s and Q^u to obtain EeiasudOi) . By expanding EeiasudOi) about the minimum energy 
orbital state for m = 0, ^i(m = 0), we obtain 

EelastiM ~ ^Kg {Oi - (m = 0))' , (51) 

where i^e = 12X\K*-K2s)/{K2sK*), and^i(m = 0) = ^cos-'[V3XA/{K2s{K*-K2s))]. By 
substituting Tc{9i) ^ T°[l + a{9i — 6i{m = 0))] and Eeiasud^i) in Eq- (^9|), and minimizing 
f{9i,m) about 9i, we obtain 

Oi{m) -9i{m = 0) = ^ \2T^am^ - r>m^l , (52) 



/(m) = 2(T - T°)m2 + -f (iT, - 2a^T;^)m^ + ©(m^). (53) 

Kg 

Therefore, we will have the second order phase transition when Kg > 2q;^T°, and the first 
order transition when Kg < 2a^T°. This result implies that when the lattice is soft (small 
Kg) or the coupling between magnetic interaction and lattice is large (large a), the transition 
becomes first-order. The estimate of the parameters for LaMnOs in Sect. VI. C predicts the 
second-order phase transition. Above equation for 6i{m) also gives the net change of 6i, 
ei{T = 0) - ^i(T > Te) = aT^/Kg. 
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IV. MODEL OF BAND STRUCTURE AND OPTICAL CONDUCTIVITY 



We use a tight-binding approximation to calculate band structure and optical conductiv- 



ity for strained films. This method is explained in detail in Ref. |T3|, and summarized in this 
section. According to band theory calculations0'El the conduction band is derived mainly 
from the Mn Cg symmetric d-orbitals and is well separated from other bands. Therefore, we 
only consider Mn Cg levels. Kinetic energy and chemical potential terms are 

Hke + Ha = -- y tfdl rf;, + H.c. - fiY dl dr,^. (54) 

i,5,a,b,a i,a,a 

Here i represents the coordinates of Mn sites, 6{= ±x,y,z) labels the nearest neighbors of 
Mn sites, a and b represent the two degenerate Mn Cg orbitals on a site, a denotes the spin 
state, and t'g' is the hopping amplitude between orbital a on site i and b on site i + 6. We 
choose {ipi >= iSz"^ — > and \ip2 >= {x"^ — y^>as the two linearly independent Cg orbitals 
on each site as before. The hopping matrix t'g' has a special form: For hopping along z 
direction, it connects only the two |3z^ — > states, thus = t"^^ = to for a=b=l, and 
zero otherwise. The hopping matrices in other bond directions are obtained by appropriate 
rotations. The Jahn- Teller coupling for uniform distortion and staggered Q2 distortion 
with wave vector Kiatuce = (tt, it, 0) or (vr, vr, n) is 

The Hund's coupling for antiferromagnetic core spin configuration with wave vector Kgpin 
((0, 0, vr) for A type AF, (vr, vr, vr) for purely AF) is 

^^Hund = -^h5'c J2 [(1 ~ e^iiKspin " i)) 4.a,T^?.«-T + + exp{iKspin ■ i)) 4.aa^?.«a ■ (56) 

i,a 

The total Hamiltonian is the sum of the terms considered so far. By diagonalizing this in k 
space, we can find the energy levels and eigenstates. 

Optical conductivity per volume due to the transitions between Mn Cg levels, a, can be 
calculated using the eigenstates and energy levels found from the above Hamiltonian. Using 
the standard linear response theory,0'§ optical conductivity is given by 




cr, 



xu _ 1 < 0\Jl^\n >< n\JpJO > 



where e is an infinitesimal and Jp is given by ^ = — ^|f^ 6 a b a^'f^ {d\aJ^i+5ba ~ H.c}j^ 

V. COMPARISON OF (7r,7r,0) AND (7r,7r,7r) TYPE DISTORTIONS 

In the purely local model considered in Ref. for which only nearest-neighbor Mn-0 
and Mn-Mn springs are considered, we have Ks = and E^^/'^ {utiit) = Ef^/^ {rniO) up to 
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order according to Eqs. p^) and (Q). However, in real situation for which farther ion- 
ion elastic energies exist, we have Kg ^ and the two distortions have different energies. 
Equations ( p2D and (|3^ ) shows that if < Kg, then the leading order term stabilizes the 
(tt, TT, 0) distortion over the (vr, vr, vr) distortion. If Kg and are very close to K*, the 

term will be smaller for the larger of Kg and K2S, opposite trend of the leading order 
term. For parameter values determined in Appendix B, we find that K2S < Kg is required to 
make (vrvrO) distortion favored. Optic phonon spectrum along k = k{7i,7i,0) and ^(Tr, 7r,7r) 
{0 < K < 1 ) would be useful to check this condition. At k = 0, these two modes have 
the same energy. As k, approaches to 1, if the (7r,7r,7r) mode has a higher energy than the 
(tt, tt, 0) mode, it would be an indication that Kg > K2g. So far, the phonon spectrum for 
LaMnOs has not been calculated. Ghosez et al^ have calculated phonon spectra for similar 
compounds, BaTiOs, PbTiOs, and PbZrOs. These results show that (tt, vr, 0) mode (M2 
point in Fig. 1 of Ref. |l^) has a higher energy than (7r,7r,7r) mode (-R12' point in Fig.l. 
of Ref. |I^) by 11 %, 7%, and 2 % respectively, contrary to our expectation for LaMnOs. 
Since the energies of these modes depend sensitively on transition metal elements, a phonon 
spectrum calculation for LaMnOs is necessary. 

We examine the possibility that the two distortions have different band energies. For 
this purpose, we use a tight binding Hamiltonian introduced in the previous section, and 
calculate the band structures and total band energies for the two distortion patterns. On 
the first Brillouin zone boundary (i.e., \kz\ = 7r/2 or + = tt planes) and on the planes 
satisfying \kx\ = \ky\, we find that the band structures for (vr, vr, 0) and (vr,vr,vr) distortions 
are identical. Between these planes, when a band is well separated from other bands, it 
has a similar band structure for the (vr, vr, 0) and (vr, vr, vr) distortions. Since the filled bands 
are well separated from the empty bands by the Jahn- Teller splitting, the results show that 
the change of the filled bands are negligible. We find that the total band energy per Mn 
ion changes less than 1 meV between the (vr, vr, 0) and (vr, vr, vr) distortions. We therefore 
conclude that the strain (lattice restoring force) effects are crucial. 

In our calculation for the strained films, we assume that the sign of the energy difference 
between the two types of ordering is not changed by apphed strains, and consider the (vr, vr, 0) 
ordering only. 



VI. RESULTS 

A. Lattice and orbital states in strained films 

In this section, we present our calculations for strained films. We use e\\ to denote the 
substrate induced additional strain, e^x — ^^xx^- We examine the strain in the range of -2 % 

o 

< e|| < 2 %. Since Oq ~ 4.03 A, the difference of the xy plane lattice parameter between 

o o 

bulk LaMnOs and substrate is between —0.08 A and 0.08 A. We represent the changes in 
Q^ui Q2S, and e^z by writing 6 in front. For example, SQ^u = Qsu — ■ Our calculation 
shows that 5ezz versus e\\ in this range is close to linear. We obtain 5ezz/^\\ ~ -1-8 for the 
parameter set obtained from Ref. 0, and ^Czz/^w ~ -1-4 for the parameter set obtained from 



Ref. 20. These ratios are about 2 times larger than the ratio for Lao.rCao.sMnOs film.E3 This 



is due to the softening of the Qsu mode in LaMnOs. 
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From Sczz versus e\\, we can also find 5Qiu and SQsu versus e\\. The results are shown in 
Fig. 2. To understand these results, we find the leading terms of SQi^/ {aoe\\), SQ2s/{cLoe\\), 
and SQsu/{cLoG\\)- They are 

^"^^ ."^J.T/".. . +0{A% (58) 



aoeii Kb + 2 {K* - 

SQ2S 3 Kb{2K*~3K2s)+A{K*-K2s) 



aoeii V2 {Kb + 2K* - 2K2s) {K* - K2s) Kl 



■\A^O{A^), (59) 



v^ .^ , o..;. . ^0{A\ (60) 



aoeii Kb ^2 {K* - K 



2s 



This shows that Q2s changes more slowly than Q3„, because the staggered distortion is 
coupled to the uniform strain only through the anharmonic term. Since K* — K2s is almost 
one order smaller than Kb-, Qiu also changes slowly compared with Q^u- Therefore, the 
main effect of the uniaxial strain is the change in the uniform tetragonal distortion without 
much change in staggered distortion or volume. Figure 2 shows that 6Q2s/{O'0^\\) can be 
either positive or negative depending on the parameter values, whereas SQiu/{aoe\\) > and 
SQ3u/{o,o^\\) < always. 

o 

We obtain 6i ^ 54.8 A for the bulk state from Eq. (^). 6i versus e\\ is shown in Fig. 



3. The change in 9i is about ±5° — 15". For tensile strains, 9i and 62 approach towards 
90°, which corresponds to — >. For compressive strains, 61 and 62 approach to 0° 
and 180°, which corresponds to |32;^ — >. This can be understood from the fact that in 
6i,62=\x'^ — y'^ > state, the x-y plane Mn-O-Mn distance tends to be farthest, and in 61, 62 
= iSz"^ — >, shortest due to the electron distribution in the xy plane. 

The change in the Jahn- Teller splitting, SAEjt, is shown in Fig. 4, which is about ± 
0.02-0.2 eV for ± 2 % strain. 



B. Magnetic property in strained films 



According to Ref. |l], J.^y = -1.66 meV, J^=1.16 meV for bulk LaMnOs. Since Oi = 54.8° 



for bulk, we obtain Jp = 7.52 meV and Jaf = 4.50 meV, which corresponds to J^y = —0.22 



meV, J" = 2.62 meV, and J* = 1.88 meV. Using these parameter values, we find Jxy and 
Jz versus 9i, which are plotted in Fig. 5 (a). As the two orbital states approach to 7r/4 
and 37r/4, i.e., {z"^ — x"^ > and — y^ >, both J^y and Jz become more ferromagnetic, 
and as they approach to and vr, i.e., {Sz"^ — > and |32^ — >, or to 7r/2 and 11/2, 
i.e., \x^ — y"^ > and — >, Jxy and Jz become less ferromagnetic. This is due to 
the orbital-state-dependent hopping between the filled and empty orbitals, which mediates 
ferromagnetic interaction. When 20° < ^1 < 70°, the magnetic ground state remains A- 
type antiferromagnetic. Outside this range, both J^y and Jz become positive, and purely 
antiferromagnetic state is the ground state. In fact. Fig. 3 shows that 2 % tensile strain can 
change 61 close to 70°, and turn the material into a purely antiferromagnetic state. 

The mean field estimates of Tc are shown in Fig. 5 (b). Tc for bulk state is about 210 
K, somewhat larger than the measured Tc = 140 K.ill It shows that ± 2 % strain changes 
Tc by about ± 50 K. Relatively large change in Tc and the possible change into purely AF 
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state are due to the strong dependence of magnetism on the Cg orbital state and the strong 
Jahn- Teller coupling between the Cg orbital state and the lattice distortion. 



C. Magnetostriction effects in bulk state 

In this subsection we present the magnetostriction effects calculated by the model in 
Section III. B. We use the two sets of parameter values in Table I and J* = 1.88 meV. The 
results obtained by numerical minimization of the Landau free energy [Eq. (^Qf)] are shown 
in Table II. The change in the JT strain, Se*, is about 0.003 ~ 0.01. 6e*/[e*{T > Tc)] is about 
-0.08 for the parameter set from Ref. ^ and -0.31 for that from Ref. [2^. These results show 



that when the effective JT modulus of Q^u niode, K* — K2s, is smaller, the magnetostriction 
effect is larger. 

We obtain a ~2, and 2a^T° ^0.1 eV for both parameter sets from Refs. ^ and We 
obtain Kq =1.2 eV for Ref. El, and Kq =0.3 eV for Ref. EO. Therefore, the transition will 



be of the second order. However, if we have a softer Kq (a third or a tenth), or a stronger 
magnetostriction coupling a (twice or three times), then we will have a first-order phase 
transition. Our numerical minimization of the free energy confirms these results. 

Recently, the orbital ordering in LaMnOs has been directly observed using a resonant 
X-ray scattering technique.il In this result, the orbital ordering versus temperature curve has 
a change of the slope at T = T/y. The sign of the change indicates that the orbital states 
change away from 3x^ — r^/3j/^— (6'i=60°) as T ^ below T^, which is consistent with our 
calculation. Recent neutron diffraction study measured sin^i versus T.il Our results predict 
that sin6'i (which corresponds to C2 in Fig. 3 in Ref. |2^) changes by 0.01 for parameters 
from Ref. |] and 0.04 for parameters from Ref. ^ between T=0 and T > T/y. However, Ref. 



22| shows negligible change in sin^i between T=0 and T > Tat, which indicates a smaller 



magnetostriction coupling a, or a larger elastic modulus Kq than the values obtained above. 



D. Band structures in strained films 

First, we summarize the band structure in bulk state A-type antiferromagnetic LaMnOs, 
which is explained in detail in Ref. |l^. Crudely speaking, the bands fall into 4 pairs, which 
may be understood by setting to = [as occurs at A; = (7r/2, 7r/2, 7r/2)]; in this case we 

have four separate energy levels on each site, which are Ei^2 = ~ \I\^\J Qls + Qlm -^3,4 = 

\[\^\JQ2s + Qlu-, ^5,6 = 2Jh'S'c - \I\^\IQ2s + Qlui ^ud E-j^^ = 2JhS'c + \[\^\JQ2s + Qlu- 

To find the three parameter values of our model Hamiltonian, i.e., t^, A, and Jh'S'c, we fit 
our band structure calculation to the LDA (local density-functional approximation) band 
calculation for the JT distorted LaMnOa in Ref. at high symmetry points in reciprocal 
space. The standard deviation is ~ 0.2 eV. The determined parameter values are to=0.622 

o 

eV, A=1.38 eV/A, and 2JhS'c=2.47 eV. The fitted band structure is shown in Fig. 1 in Ref. 
13|. 

When the strain does not change A-type antiferromagnetic core spin configuration, the 
main effect of the strain is the change in the band width. The results are shown in Fig. 6. 
Solid lines are for the compressive strain, and dashed lines are for the tensile strain. Bulk 
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band structure can be approximately obtained by taking the average of the two band struc- 
tures. At (7r/2, 7r/2, '7r/2), where the effective hopping vanishes, the energy level change is 
6AEjT obtained before. From this point, dispersions along {n/2, 7i/2, 0) and along (tt, 0, 7r/2) 
represent the hoppings in z direction and in xy plane, respectively. Between (7r/2, 7r/2, 7r/2) 
and (7r/2, 7r/2, 0), the widths of the lower JT bands, Ei^2 and -^5,6, are increased for com- 
pressive strains, whereas the widths of the upper JT bands, £'3,4 and -^7,3, are decreased. 
This is related to the changes in 9i and 62 due to the strains: As 9i and 62 approach to 
and TT, the lower JT level approaches to 3z^ — state which has a large hopping along z 
direction, whereas the upper JT level approaches to x"^ — which has no hopping along z 
direction. But between {n/2, 7r/2, 7r/2) and (tt, 0, tt/2), the dispersion does not change much, 
indicating that the average hopping is not changed much in xy plane due to the alternating 
orbital pattern in xy plane. 

As pointed out in the previous section, 2 % tensile strain may induce purely antiferromag- 
netic ground state. Due to the strong Hund's coupling, this results in substantial changes 
in band structure and optical conductivity as shown below and in the next section, respec- 
tively. Band structure for the purely AF state can be obtained by using Kspm = (tt, tt, tt) in 
the model described in Sect. IV. The results are shown in Fig. 7 for the same lattice distor- 
tions and parameters used for the band structure shown as dotted lines in Fig. 6. Between 
(7r/2, 7r/2, vr/2) and (7r/2, 7r/2, 0), the two band structures are identical since it involves only 
z direction hopping. Between (7r/2, 7r/2, 7r/2) and (7r,0,7r/2), where only xy directional 
hopping is involved, for A-type AF state the different JT levels repel each other, while for 
purely AF state the different Hund's levels repel each other. This represents different mix- 
ings for different spin configurations: the mixing for A-type AF state is mainly between the 
different JT states, while the mixing for purely AF state is mainly between the different 
Hund's states, the band structure between (7r,0,0) — (0,0,0) — (7r/2, 7r/2, 7r/2) shows the 
same trend, which results in a small indirect band gap for the purely AF state. However, 
due to the on-site Coulomb repulsion of about 1-2 eV neglected in the above calculation!!^ 
it is unlikely to have the insulator-to-semimetal transition in this material. 

E. Optical conductivities in strained films 

From the band structure, we have calculated optical conductivities for strained films in 
A-type AF ground state. Results are shown in Fig. 8. a^x or ayy (solid lines) shows a 
relatively small changes by strains. The spectral weight of the Hund's peak in (dotted 
lines) at around 2.5 eV is increased (decreased) about 20 % as we apply 2 % compressive 
(tensile) strain. This seems consistent with the change of the average hopping along z 
direction by strains mentioned in the previous section. 

Optical conductivities for the purely AF state are calculated in the same way. Figure 9 
shows the results for a^x (solid line) and a^z (dotted line). First, the sharp Hund's peak in 

is disappeared. This can be understood from the band structure, particularly between 
{71/2, 7r/2, Tc/2) and (tt, 0, 7r/2). The sharp Hund's peak in (Xzz for A-type AF state originates 
from the two parallel bands split by the Hund's coupling, 2Jh5'c. However, this structure 
disappears when we have purely AF core spin as seen in Fig. 7. Comparison of Fig. 9 (a) (2 
% strain) and (b) (bulk) shows that as orbital state is changed toward x"^ — y"^ by the tensile 



15 



strain, which has zero hopping along z-direction, the spectral weight of the Hund's peak in 
Gzz decreases, as observed for the A-type AF state in Fig. 8. a^x shows prominent peaks at 
the Hund's splitting (2.5-3.5 eV) and at the Hund-plus-JT splitting (4-5 eV) due to purely 
AF spin state in contrast to the A-type AF state (see Fig. 8). The JT peaks at around 1 
eV in both Oxx and are due to the strong hybridization between major and minor spin 
states. If we increase JhSc and reduce this hybridization, the JT peaks decrease, as can be 
seen by comparing Figs. 9(a) and 9(c). Above results show that changes in xy plane spin 
configuration make differences in a^z due to the changes in hybridization. 

For the calculations so far, we have assumed zero on-site Coulomb repulsion \J . However, 
as explained in Ref. 0, in this material there exists f/ ^ 1.6 eV. Therefore actual peak 
positions will be higher by ~ 1.6 eV and the spectral weight of each peak will be reduced 
inversely proportional to the peak energy. 



In summary, we developed a model of elastic energy for LaMnOa and solved for uniaxial 
strains in thin films. We found that ±2 % strain can change the uniform tetragonal strain and 
tg orbital states without much change in the staggered distortion or volume. We found that 
2 % tensile strain can change the magnetic ground state into purely antiferromagnetic state, 
inducing dramatic changes in band structure and optical conductivity. Magnetostriction 
effect at Tn in bulk state is found to be large. We examined the possibility that the lattice 
energy will favor (vrvrO) ordering over (tttttt) ordering. We also noted "Qsm mode softening" 
in LaMnOs. The results presented in this paper for LaMnOs show the strong coupling 
between lattice, Mn Cg orbital state, and exchange interaction, which lies at the root of the 
novel properties of doped manganese perovskites. 

This work is supported by NSF-DMR-9705182 and the University of Maryland MRSEC. 

APPENDIX A: GENERAL EXPRESSION OF ELASTIC ENERGY FOR 

STAGGERED DISTORTION 

In this Appendix, we show how we can get the general expression of the elastic energy 
due to distortions with wave vector k in perovskite structure. We again consider three 
dimensional displacement of Mn ion, 5, and the displacements of oxygen ions along Mn-Mn 
axis, M^, M^, and as defined in the text and shown in Fig. 1. We consider the displacements 



with wave vector fc, 5^ = 5re*'^'* and wi'^'^ = w^'^'^e*'^'*. When we define o s(^) = 5?'^'^ and 



In special cases, symmetry arguments make certain terms vanish or equal. First, when 
kx = vr, mirror symmetry operation about x = plane changes kx to —kx- Since kx = tt 
and k'^ = —n are equivalent, the symmetry operation changes k back to k. Since this 



VII. CONCLUSION 



4,5,6(^) = U 




(Al) 



where Dij{k) = Dji{k). 
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— * 

operation changes 5| to —55, odd order terms of (5| should vanish, and therefore, Di2{k) = 
Di^{k) = Di4(k) = Di^{k) = DiQ^k) = 0. Second, when k^ = 0, mirror symmetry operation 
about X = plane changes k back to k, 6f to —Sf, and uf to —uf- Therefore, odd order 
terms of 5| and vanish. Therefore, Dui^k) = Di^{k) = Di^{k) = DiQ^k) = and 
Di2{k) = D4^{k) = Di^{k) = DiQ{k) = 0. Third, when k^ = ^ky, mirror operation 
interchanging x axis and ±y axis changes k back to k, and 6^ to ±5y, to ±.Uy. Therefore, 
Diiik) = D22{k\ D,^{k) = D,,{k), Du{k) = D^^ik), D^^ik) = D^^ik), Dij{k) = ±D2j{k), 
Dij{k) = ±D5j{k), where j = 3, 6 

If we apply these rules to k = {it, it, 0), and (vr, vr, tt) distortions, we obtain the following 
expressions. 

Dn{k){6f + 5f) + D^,{k)5f + D,,{k){uf + uf) 
+D,S)uf + 2D^,{k)ulul + 2D;,{k)5lul. (A2) 

D,,{k){5f + 5f + 5f) + D,,{k){uf + uf + uf) 

+2D,,{k){ulul + ulul + ulul). (A3) 

K2S1 K^s, and Kg in the text are defined as K2S = [-D44(7r7r0) — D45(7r7r0)]/2, K^g = 
3[D44(7r7rO) + D45(7r7r0)]/2, and Ks = [Z^44(7r7r7r) — 04^^(717171)] /2. Therefore, the condition 
for K^s > K2S is D44(7r7r0) > — 2D45(7r7rO). We expect this condition is well-satisfied for the 
following reasons: First, the stability of lattice implies Z?44(7r7r0) > 0. Second, the Coulomb 
repulsion between oxygen will favor similar 0-0 distances, which implies D45(7r'7r0) > 0. 



£;,[fe = (7r,7r,0)] 



— * 

Es{k = (tt, 7r,7r)] 

Nmu 



APPENDIX B: PARAMETERS FOR THE MODEL OF ELASTIC ENERGY 

Optical absorption experiment for LaMnOs shows Mn-0 bond stretching mode peak at 

70.3 meV.ii From this we can find the effective Mn-0 spring constant -ft'i=7.36 eV/A (Ref. 
O). If we assume that Ki is the main contribution to K2S1 we obtain K2S ~ Ki/2 = 3.68 

o2 



eV/A . In Ref. |T3|, we used this value of Ki to estimate the Jahn- Teller coupling constant 
A, the value of which is consistent with LDA band calculation.lli 

For Lao ssSro lyMnOs, cn and C12 have been measured in an ultrasound experiment.S 

We believe the bulk modulus cb = (cn + 2ci2)/3 =143 GPa(is:ij=10.8 eV/A ) of 
Lao.ssSro.irMnOs at 200 K (orthorhombic structure) can be used as an approximate value 
for the bulk modulus of LaMnOs. The value of cb at 310 K (rhombohedral phase) is about 
176 GPa, which gives a rough estimate of uncertainty in cb of about 20 %. 

Since Lao.ssSro.iyMnOs is in the doping region where a structural change happens,^ and 
c* is sensitive to the structural transition unlike bulk modulus, we do not believe the mea- 
sured c* of this material is a good estimate of c* for LaMnOs. Indeed, c* for Lao.ssSro.irMnOs 
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is much smaller than those for other typical perovskite oxides.E^ For example, the LDA calcu- 
lation in Ref. H predicted cb = 199 GPa, c* = 142 GPa, and c*/cb = 0.71 for SrTiOg, which 
is close to the measurements for SrTiOa, cb =179 GPa, c*=115 GPa, and c* /cb = 0.64. The 
same LDA calculation results showed that other perovskite oxides, such as BaTiOa, CaTiOa, 
KNbOs, NaNbOs, PbTiOg, PbZrOg, BaZrOsHowever, for Lao.ssSro.irMnOs, cb=U3 GPa, 
c*=48 GPa , and c*/cb=OM at 200 K and cb=176 GPa, c*=35 GPa , and c*/cb=0.20 at 
300 K. Therefore c*/cb ratio is less than half of the values for typical perovskites. 

Therefore, we instead use the C12/C11 ratio measured in Lao.rCao.sMnOs thin film (Ref. 
^ and Lao,6Sro.4Mn03 thin film (Ref. since their doping ranges are relatively far from 
the structural phase transition doping ratios. The results are Ci2/cii=0.312 for Ref. and 



ci2/cii=0.374 for Ref. ||. Using 

c* 3 1 - cu/cu 



Cb 2 1 + 2ci2/cii' 



(Bl) 



we obtain c* and K* shown in Table I. 

o 

We obtain Q2S QsL for bulk LaMnOa from crystallography data: (52s=0.398 A, 

o 

(^3^ =-0.142 A (Ref. |[). Therefore, the three unknown quantities of the model, Qll, A, and 
A, are determined from the equilibrium condition: 



dE 



dQiu 



dE 



0, 



eq ^^3« 



(B2) 



eq 



Obtained parameter values are shown in Table I. The values of A are small enough to 
justify our approximation of anharmonic terms. For example, the largest dropped term, 
AQlJ{?,^), is 0.5-2 % of K*Q\J2 for parameters in Table I. We obtain Q^^ = 0.024 

o o 

A and 0.005 A for the parameters from Refs. |^ and respectively, which shows that the 
average bond length does not change much. Therefore, it is reasonable to approximate by 

o 

the average Mn-Mn distance observed in bulk LaMnOs, 4.03 A. The A values are close to the 

o 

independent estimate A=1.38 eV/A obtained from band structure fitting. Negative sign of A 
is consistent with thermal expansion. It is noteworthy that the size of A is largely different 
for the two parameter sets, which is the consequence of the Qs^ mode softening: a small 

O 2 o 2 

difference in K* gives a quite large difference in K* — K^s (0.89 eV/A and 0.19 eV/A for 
the parameters from Refs. ^ and |2y, respectively), which results in a large difference in the 
estimate of A. 
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FIGURES 



FIG. 1. Degrees of freedom considered in our model, 6^ and ui'^'^ . 

2 2 

FIG. 2. 6Qiu, SQ2s, SQ3u versus e\\ for (a) Kb=10.8 eV/A , K*=3.87 eV/A , A=1.25 eV/A, 
ir2s=3.68 eV/A , ^=-0.329 eV/A , which are the parameter values from Ref. 2, and (b) ii'B=10.8 

eV/A , if* =4.57 eV/A , A=1.13 eV/A, K2s=S.68 eV/A , ^=-1.65 eV/A , which are parameter 
values from Ref. 20, shown in Table I. 

FIG. 3. 01 versus e||: The solid line is for the parameter values from Ref. 2, and the dashed 
line is for the parameter values from Ref. 20, shown in Table I. 

FIG. 4. SEjT versus e||. Solid line is for the parameter values from Ref. 2, and dashed line is 
for the parameter values from Ref. 20, shown in Table I. 

FIG. 5. (a) Magnetic coupling constant Jxy (solid line) and Jz (dotted line) versus orbital 
state 01. (b) Mean field estimation of Tc versus orbital state 0i. the ground state is A type 
antiferromagnetic between 20° and 70°, and purely antiferromagnetic outside this range. 

FIG. 6. Band structures for 2 % compressive strain (solid lines) and 2 % tensile strain (dotted 
lines) for the parameter values from Ref. 2, shown in Tabic I. Tight binding model parameter 
values are to=0.622 eV, 2JhSc=2.4:7 eV, and A=1.13 eV. A-type antiferromagnetic ground state 
is assumed. 

FIG. 7. Band structure for purely antiferromagnetic core spin state for 2 % strain for the 
parameter values from Ref. 2 shown in Table I, to be compared with the dotted lines in Fig. 6. 
The tight binding parameter values are to=0.622 eV, 2JhSc=2.47 eV, and A=1.13 eV. 

FIG. 8. Optical conductivities, (solid lines) and a^z (dotted lines), for 2 % compressive 
strain (top panels), bulk (middle panels), and 2 % tensile strain (bottom panels) for the parameter 
values from Ref. 2 (left panels) and for the parameter values from Ref. 20 (right panels) shown in 
Table I. A-type antiferromagnetic core spin state is assumed. 

FIG. 9. Optical conductivities, axx (solid lines) and azz (dotted lines), for purely antifer- 
romagnetic core spin configuration with (a) 2 % tensile strain, (b) bulk, and (c) a larger Hund 
coupling and 2 % tensile strain. Lattice distortions in Fig. 1 (a), to=0.622 eV, 2JhSc=2A7 eV, 
and A=1.13 eV are used for the calculation. 
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TABLES 



TABLE L Obtained values of parameters: cb and Kb are from Ref. 24, and K2s is from the 
optic phonon mode frequency. We use two experimental results in Refs. 2 and 20 to obtain the 
two sets of parameters shown here. Details about how these numbers are obtained are shown in 
Appendix B. 



Source 


CB(GPa) c*(GPa) 


2 2 2 3 

KB{eY/A ) K*{eY/A ) K2s(eV/A ) A(eV/A) A(eV/A ) K* -K2s{eY/A 


Ref. 1 
Ref. 20 


143 90.8 
143 76.8 


10.8 4.57 3.68 1.13 -1.65 0.89 
10.8 3.87 3.68 1.13 -0.329 0.19 



TABLE II. Magnetostriction effect (numerical calculation) 



T 


parameters 


Qln(A) 


Q2s{h 




Oi 


e* 


T=0 


from Ref. | 


0.0246 


0.401 


-0.131 


53.6° 


-0.0328 




from Ref. ^ 


0.0051 


0.411 


-0.099 


51.5° 


-0.0248 




from Ref. | 


0.0244 


0.399 


-0.142 


54.8° 


-0.0355 




from Ref. ^ 


0.0049 


0.399 


-0.143 


54.8° 


-0.0358 


Differences between 


from Ref. | 


2.5 xlO"'^ 


0.0024 


0.011 


-1.2° 


0.0028 


T = and T > Tc 


from Ref. ^ 


2.6 xlO""^ 


0.012 


0.044 


-3.4° 


0.011 
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Fig. 1, Ahn and Millis 




Fig. 2, Ahn and Millis 




Fig. 3, Ahn and Millis 




Fig. 4, Ahn and Millis 




Fig. 5, Ahn and Millis 
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Fig. 6, Ahn and Millis 
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Fig. 7, Ahn and Millis 
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Fig. 8, Ahn and Millis 
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